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Abstract We inspect the light and radial-velocity 
curves of the eclipsing binary AV Del. Using the most 
new version of Wilson - Devinney (1971) code (for more 
information please see Wilson - Wan Hamme, 2003), 
the absolute elements, fundamental orbital and phys- 
ical parameters of the system are determined. Using 
these parameters the configuration of the system are 
presented. Then, using the new SHELLSPEC code, 
we have studied and presented an accretion disc model 
for the system. We found AV Del is a semi-detached 
system in which an accretion disc is surrounded the 
primary star. Combining the radial-velocity and the 
light curve analysis, we have derived accurate abso- 
lute masses for the components of Mi=1.449 M Q and 
M 2 =0.687 M Q and radii of i? x =2.61 R Q and i? 2 =4.21 
Rq, as well as effective temperature of T 2 = 4281±15 
K for the secondary. Also, we have derived a tempera- 
ture of T=5700 K for the disc. Finally, our results are 
compared with those of previous authors. 
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1 Introduction 

AV Del (a 2 ooo = 20 ,l 45 m 31 s .47, 62000 = 11°10'26".4, 
V— 11.8) is a rather neglected variable star discovered 
by Hoffmeister (1935) and given the original designa- 
tion 184.1930. It is a suspected member of the cool 
Algols (Popper, 1996) and was reported by the au- 
thor as a double-lined system with a period of 3.85 
days and a mean spectral type of G5, although other 
sources have classified it as F8 (e.g., Halbedel 1984). 
Mader et al. (2005) presented new spectroscopic and 
BVRI photometric observations of AV Del. They deter- 
mined the radial-velocities for both components, using 
the two-dimensional cross-correlation technique TOD- 
COR (Zucker & Mazeh 1994). They showed the system 
to be most likely semi-detached, with the less massive 
and cooler star filling its Roche lobe. 

As we will see in section 2, there are evidences for 
a mass transfer in AV Del. So we thought it is useful 
to study the existence of mass transfer and the nature 
of the system using the new SHELLSPEC code which 
is a strong program to analyze those stars have accre- 
tion disc and circumstellar material. First, we present 
the analysis of the system using the Wilson - Devinney 
(1971) WD code in sections 3 and 4. We present here a 
new light and radial- velocity curves analysis, and then 
we determine the physical and absolute parameters of 
the system. The absolute parameters are given in sec- 
tion 5. Analysis of the system using the SHELLSPEC 
code is explained in section 6. 

2 Assumptions 

In this paper, we consider B, V, R and / filters obser- 
vations of Mader et al. (2005). The latest version of 
WD code is applied to solve the light and radial- velocity 
curves. After obtaining the absolute parameters of the 
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components, we used the SHELLSPEC code to derive 
the parameters of the disc. Throughout this paper, the 
subscripts 1 and 2 refer to the primary (hotter) and the 
secondary (cooler) components, respectively. 

During using WD code, for both components of 
the system, we used bolometric linear, logarithmic and 
square root law for the limb darkening and the best re- 
sult was obtained for bolometric logarithmic limb dark- 
ening law of Klinlesmith & Sobicski (1970) with the 
following form: 

I = I (l — x + xcos9 — y cos 8 ln(cos 6)) , (1) 

where the limb darkening coefficients x and y for both 
components were fixed to their theoretical values, inter- 
polated using Van Hamme's (1993) formula which are 
tabulated in Table 1. 

The gravity darkening exponent from Lucy (1967) 
and the bolometric albedos from Rucinski (2001) were 
chosen for convective envelopes (g=0.32, A — 0.5), 
which are in good agreement with the final surface 
temperature. In order to reduce the number of the 
free parameters, these parameters were kept constant 
during all iterations. Also, it is assumed that this bi- 
nary system has zero orbital eccentricity (e = 0.0) and 
that its rotational and orbital spins are synchronous 
(F\ = F 2 = 1.0). Also, black body radiation models 
are employed, and we assume that there is no third 
light (?3 = 0.0) for the system. 

Mader et al. (2005) calculated a linear ephemeris: 

Min.I(HJD) = 2,450,714.34779(39) 

+3.8534528(35)^, (2) 

where E is the number of cycles counted from the epoch 
of reference and the uncertainties are given in parenthe- 
ses in units of the last significant digit. A period study 
of AV Del was carried out by Qian (2002) on the basis 
of fewer times of eclipse. A small period increase of 
about 3.15 ± 0.19 x 10~ 6 days yr^ 1 was reported by 
Qian (2002), while the other quadratic fit derived by 
Mader et al. (2005) gives: 

Min.I(HJD) = 2,450,714.34771(34) 
+3.8534620(42)^ + 0.50(16) x 10~ 8 £; 2 , (3) 

which corresponds to a period change of 0.95 ± 0.30 x 
10~ 6 days yr^ 1 . These period changes may be related 
to the mass transfer presumably taking place in the sys- 
tem. The main goal of this paper is that show whether 
these period changes are related to the mass transfer or 



not? Also, we want to find if there is a real mass trans- 
fer in the system, what kind of material is corresponded 
it? 



3 Photometric Solutions 

We used the photometric data which reported by Mader 
et al (2005) for our solutions. We applied WD code for 
photometric solutions. Solutions are found for the in- 
clination angle (i), the mass ratio (q = M 2 /Mi), the 
modified gravitational potentials (f2, in the WD usage), 
the relative monochromatic luminosity of the primary 
(Li) and the mean temperature of secondary (T 2 ). The 
temperature of the primary was held fixed at the spec- 
troscopic value of Ti=6000K. The monochromatic lu- 
minosity of the secondary (L 2 ) was computed by the 
program directly from the temperatures of the primary 
and secondary, the luminosity of the primary, the radi- 
ation laws, and the geometry of the system. Other WD 
code parameters kept fix along with solutions. 

We initially obtained solutions in a mode appropri- 
ate for semi-detached systems, following the expecta- 
tion from the characteristics of AV Del. Specifically, 
we used "mode 5" for the secondary which is filling its 
Roche lobe. We also tried to use the detached configu- 
ration to complete our work. But, similar to Mader et 
al.'s (2005) results, our best results determined for the 
semi-detached configuration. 

Initially, the light curve program (LC) of WD code 
was implemented to employ an initial set of values for 
the next step of solutions which was using the differen- 
tial corrections main program (DC) of WD code. After 
several runs of DC program, we obtained the final set 
of values. The results of photometric solutions with 
the final elements are given in Table 2. The theoretical 
light curves computed with these results (continuous 
line) for each filter are shown in Fig. 1. 

4 Spectroscopic Solutions 

Using the radial- velocity data from Mader et al. (2005), 
we obtained the orbital elements of the system. Similar 
to photometric solution, we used WD code for spec- 
troscopic solution. Solutions are found for the orbital 
semi-major axis a and the radial-velocity of the binary 
system center of mass V 1 . Also, the ephemeris adopted 
is that in equation (2). The results of spectroscopic so- 
lutions are given in Table 3. Using the final elements of 
the object, the theoretical radial-velocity curves (con- 
tinuous lines) are shown in Fig. 2. 
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Table 1 The limb darkening coefficients for AV Del 



Parameters xi(bol) X2(bol) y 1 (bol) y2(bol) x\ y\ X2 2/2 



Filter B 

Filter V 

Filter R 

Filter I 



.639 0.632 

.639 0.632 

.639 0.632 

.639 0.632 



0.225 0.150 

0.225 0.150 

0.225 0.150 

0.225 0.150 



0.739 0.253 

0.823 0.165 

0.647 0.272 

0.554 0.265 



0.813 -0.038 

0.842 -0.232 

0.750 0.092 

0.652 0.150 





Fig. 1 The observed and theoretical light curves of AV Del, while the theoretical light curves are drawn on the basis of the 
values derived via WD code. Solid circles show the observed data and the theoretical light curves are shown by continuous 
lines. Top right and left diagrams show / and R filters, respectively. Also, bottom right and left diagrams show V and B 
filters, respectively. 
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Table 2 The photometric parameters of AV Del 



Parameters This work This work This work This work This work Mader 

B Filter V Filter R Filter I Filter BVRI Filters et al. (2005) 



i 


81.395 


80.543 


81.333 


81.530 


81.162 


81.42 




± 0.221 


± 0.152 


± 0.127 


± 0.167 


± 0.091 


± 0.10 


<1 


0.4973 


0.4861 


0.4916 


0.4609 


0.4907 


0.4874 




± 0.0131 


± 0.0125 


± 0.0108 


± 0.0159 


± 0.0068 


± 0.0046 


Ti 


6000 


6000 


6000 


6000 


6000 


6000 




fixed 


fixed 


fixed 


fixed 


fixed 




T 2 


4245 


4239 


4260 


4327 


4258 


4264 




± 21 


± 16 


± 15 


± 22 


± 9 


± 13 


Oi 


5.7218 


5.6405 


5.8179 


5.5627 


5.8790 


5.54 




± 0.1275 


± 0.1119 


± 0.1191 


± 0.1405 


± 0.0482 


± 0.11 


2 


2.8707 


2.8490 


2.8597 


2.7998 


2.8580 


2.852 














± 0.015 


Li/(Li + L 2 ) 

i_ / \ i_ \ • a — m a j 


0.7697 


0.6969 


0.5985 


0.5452 


0.6363 




L 2 /(Li + Lo) 


0.2303 


0.3031 


0.4015 


0.4548 


0.3637 




log gi(CGS) 


3.78 


3.77 


3.80 


3.77 


3.77 


3.759 


log g 2 (CGS) 


3.03 


3.03 


3.03 


3.03 


3.03 


3.032 


ri (pole) 


0.1911 


0.1937 


0.1875 


0.1957 


0.1853 


0.1974 




± 0.0047 


± 0.0042 


± 0.0042 


± 0.0054 


± 0.0016 




T\ (point) 


0.1933 


0.1960 


0.1895 


0.1980 


0.1872 


0.1999 




± 0.0049 


± 0.0044 


± 0.0044 


± 0.0057 


± 0.0017 




ri(side) 


0.1921 


0.1947 


0.1884 


0.1968 


0.1862 


0.1985 




± 0.0048 


± 0.0043 


± 0.0043 


± 0.0055 


± 0.0017 




ri(back) 


0.1931 


0.1957 


0.1892 


0.1977 


0.1870 


0.1996 




± 0.0049 


± 0.0044 


± 0.0044 


± 0.0056 


± 0.0017 




r 2 {pole) 


0.2993 


0.2976 


0.2984 


0.2935 


0.2983 


0.2978 




± 0.0021 


± 0.0020 


± 0.0017 


± 0.0027 


± 0.0011 




r 2 (point) 


0.4287 


0.4264 


0.4275 


0.4211 


0.4274 


0.4267 




± 0.0085 


± 0.0083 


± 0.0071 


± 0.0112 


± 0.0045 




r 2 (side) 


0.3124 


0.3105 


0.3115 


0.3061 


0.3113 


0.3107 




± 0.0022 


± 0.0021 


± 0.0018 


± 0.0029 


± 0.0012 




r 2 (back) 


0.3449 


0.3430 


0.3440 


0.3387 


0.3438 


0.3433 




± 0.0022 


± 0.0021 


± 0.0018 


± 0.0028 


± 0.0011 






2.8707 


2.8490 


2.8597 


2.7998 


2.8580 




^out 


2.5735 


2.5579 


2.5656 


2.5224 


2.5643 




E uj(o-c) 2 


9 x 10~ 4 


5 x 10~ 4 


3 x 10" 5 


4 x 10" 5 


4 x 10~ 4 
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Table 3 The spectroscopic parameters of AV Del 



Parameters This Work Mader 

et al. (2005) 



V^Km/s) -64.141 -63.14 

± 0.181 ± 0.24 

a(i? Q ) 13.3192 13.34 

± 0.0347 ± 0.09 




Fig. 2 The observed and theoretical radial-velocity curves of AV Del, while the theoretical radial- velocity curves are drawn 
on the basis of the values derived via WD code. Solid circles show the observed data and the theoretical radial-velocity 
curves are shown by continuous lines. 
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5 Absolute Elements 



6 SHELLSPEC Solutions 



Using the obtained results of the light and radial- 
velocity curves, we calculated the absolute parameters 
of the system. The results are listed in Table 4. We de- 
termined the absolute dimensions using the following 
formulae: 



Mi/Mq — fi(Mi, M 2 ,i)(l + q) / sin i, 



f 1 (M 1 ,M 2 ,i) = (M 2 smif/(M 1 +M 2 f 



M 2 /M Q =q(M 1 /M Q ), 



(4) 



(5) 



(6) 



Ri, 2 /R Q = 4.207(M l!2 /M Q (l + q) P 2 (day s)) 1/3 r (side), 

(7) 



L h2 /L Q = (R h2 /R Q ) 2 (T h2 /T Q )\ 



(8) 



P1.2IPQ = (0.01344Mi )2 )/[(Mi + M 2 )P 2 rl 2 l (9) 



(M w )i, 2 = 42.36-lOlogT 1 , 2 -51og(i? 1 ,2AR ). (10) 

The fillout for the components can be calculated 
from the following formula: 



fillouti, 2 = x 100. 

"1,2 



(11) 



Also from the following formulae we calculate K~i >2 
and a 1,2 sin i for each of the systems: 



if 1,2 = 



27rai, 2 sini 



. . , GP 2 h{M u M 2 ,i) ^ 
aismi = [ — 2 ]s : 



a 2 sin i = 



a\ sini 



(12) 



(13) 



(14) 



The results are tabulated in Table 5. On the basis 
of these values, we have drawn the configuration of the 
components using the Binary Maker 2.0 (Bradstreet, 
1993) software, which are shown in Fig. 3. 



Assumption the third body light in the WD code solu- 
tions (^3 ^ 0.0) was not terminated to the expedient re- 
sult, but it should be mentioned that the WD code does 
not consider the complete formulae about nor accretion 
disc, nor another environmental body. So, it is neces- 
sary that we use a more complete code (or program) to 
solve the light curve of the binary systems that may be 
surrounded by a third body. The SHELLSPEC code is 
designed to solve simple radiative transfer along the line 
of sight in three-dimensional moving media. The scat- 
tered light from a central object is taken into account 
assuming an optically thin environment. Output inten- 
sities are then integrated through the two-dimensional 
projection surface of a three-dimensional object. The 
assumptions of code include LTE and optical known 
state quantities and velocity fields in three dimensions 
(Budaj et al. 2004). Now, using the results obtained 
in sections 3, 4 and 5 we seek the parameters of a third 
body in the system with the SHELLSPEC code. In 
this code both objects, the primary and the compan- 
ion (secondary), may have shapes according to Roche 
model for detached, semi-detached or contact systems. 
So, according to previous sections, we assume a Roche 
model for semi-detached configuration that the primary 
is a nontransparent uniformly rotating sphere and black 
body radiation models are employed to solutions using 
the SHELLSPEC code, too. 

The limb and the gravity darkening coefficients used 
in this code are those explained in section 2. But the 
reflection effect is not included in calculation. 

Since there are too much free parameters in the 
SHELLSPEC code (because there are several environ- 
mental body), in this stage, we kept the parameters 
of the primary and the secondary (those obtained in 
sections 3, 4 and 5) fix and studied the effect of the pa- 
rameters related to the accretion disc, spot, jet, stream 
and shell on the light curve to obtain the best fit to the 
observed data. 

Keep fixing the free parameters of the primary and 
secondary stars on the values obtained on the basis of 
the WD code solutions will not eventuate a bad result, 
because the WD code, noticeably is a perfect code to de- 
rive the parameters of the stars those are in the binary 
systems and thus we will not centralize our calculation 
on the free parameters of the primary and secondary 
and will attend to the other objects. 

Beside, we only consider R filter observation of 
Mader et al. (2005) for our solutions, since as it is 
shown in Table 2, the accuracy of the results for this 
filter are more than others. 

Solutions were determined for dine (angle between 
rotation axis of the model and the line of sight), dd 
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Table 4 The absolute elements of the binary system AV Del 



Parameters This Work This Work This Work This Work This Work Mader 

Filter B Filter V Filter R Filter I Filter BVRI et al. (2005) 



M 1 /M Q 


1.426 


1.437 


1.432 


1.462 


1.433 


1.453 


M 2 /M Q 


0.709 


0.699 


0.704 


0.674 


0.703 


0.705 




6.340 


6.489 


6.094 


6.640 


5.950 




L 2 /L Q 


4.399 


4.333 


4.440 


4.572 


4.432 




Pi/Pe 


0.085 


0.083 


0.091 


0.081 


0.094 




P2/P& 


0.010 


0.010 


0.010 


0.010 


0.010 




Ri/Rq 


2.56 


2.59 


2.51 


2.62 


2.48 


2.632 


R2/R0 


4.26 


4.24 


4.25 


4.18 


4.25 


4.233 


Mi (bol) 


2.58 


2.56 


2.63 


2.53 


2.65 


2.47 


M 2 (bol) 


2.98 


3.00 


2.97 


2.94 


2.97 


2.90 


h{M l ,M 2 ,i)/M Q 


0.076 


0.072 


0.074 


0.065 


0.074 





Table 5 Other elements of the binary system AV Del 



Parameters This Work This Work This Work This Work This Work Mader 

Filter B Filter V Filter R Filter I Filter BVRI et al. (2005) 



a\ smi/ Rq 


4.371 


4.292 


4.332 


4.149 


4.332 


a 2 sin i / Rq 


8.789 


8.829 


8.812 


9.002 


8.828 


K^Km/s) 


57.434 


56.396 


56.922 


54.517 


56.922 


K 2 (Km/s) 


115.486 


116.012 


115.788 


118.285 


115.999 


fillout 1(%) 


50.171 


50.510 


49.153 


50.332 


48.614 


fillout 2(%) 


100.000 


100.000 


100.000 


100.000 


100.000 
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Fig. 3 The Configuration of the binary system AV Del on the basis of the values derived via WD code and drawn using the 
Binary Maker 2.0 (Bradstreet, 1993) software. From top to bottom, first, second and third figures show the configuration 
of the system in Phases 0.00, 0.50 and 0.75, respectively. And the last figure shows the configuration of the system with its 
Roche lobe. 
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Table 6 The parameters of the accretion disc of AV Del 



Parameters 


This work 


adisc (Rq) 


1.2 


rindc (Rq) 


3.5 


routdc (Rq) 


8.0 


tempdc (K) 


5700 


densdc (gr/cm 3 ) 


33 x 10~ 15 


anedc (cm -3 ) 


21 x 10 9 


edendc 


-1 


vtrbdc (Km/s) 


90 



(distance to the shell from the Earth), rstar (radius of 
the primary), Tstar (effective temperature of the pri- 
mary), emstar (mass of the primary), dlst (limb dark- 
ening coefficient of the primary), dgst (gravity dark- 
ening coefficient of the primary), rep (radius of the 
secondary), tempep (effective temperature of the sec- 
ondary), dlcp (limb darkening coefficient of the sec- 
ondary), dgcp (gravity darkening coefficient of the sec- 
ondary), qq (mass ratio, secondary to primary), rsp 
(radius of the spherical spot), tempsp (effective tem- 
perature of the spot) , adisc (half of the thickness of the 
disc cylinder), rindc (inner radius of the disc), routdc 
(outer radius of the disc), tempdc (effective tempera- 
ture of the disc), densdc (mean density of the disc), 
anedc (the electron number density of the disc), edendc 
(the exponent of the power-law behavior of the densities 
p ~ r v in the disc ) and Vtrb (microturbulence) . Where 
dine, rstar, emstar, dlst, dgst, rep, tempep, dlcp, dgcp 
and qq adopted from the values determined in previous 
sections. Also, Tstar and dd kept fix with values in 
Mader et al. (2005). 

First, we ran the code assuming no third body, but 
we could not obtain good fit on the observed data. Then 
we tried assuming only disc as a third body. On the 
basis of this configuration agreement between the ob- 
served and the theoretical light curve was very good, 
but we tried to exam the effect of other bodies (spot, 
jet, stream and shell) on the light curve. Finally, the 
best fit obtained while we considered only an accretion 
disc. As we see in Fig. 4, when the disc is assumed, the 
agreement between the observed data and the theoreti- 
cal light curve is quite good. The parameters of the disc 
are given in Table 6. The configuration of the system 
while its primary star is surrounded by the accretion 
disc on the basis of the values of Table 6 are shown in 
Fig. 5. 



7 Conclusion 

Comparing the new photometric-spectroscopic solu- 
tions of the systems with those of given in literature, 
the present authors can get following conclusions: 

In our study, we may conclude that AV Del is an 
eclipsing binary system which we have employed a new 
tools called SHELLSPEC code to solve its light curve. 
We have presented a new configuration for the system 
that consists an accretion disc surrounding the primary 
star. As it is given in Table 6, the accretion disc has 
the shape of a cylinder. Its vertical half-thickness (1.2 
Rq) is comparable to the radius of the primary star. 
The distance between the primary and the secondary 
centers is about 13.3 Rq and the outer radius of the disc 
is 8 Rq. So, since the radius of the primary and the 
secondary are about 2.6 Rq and 4.2 Rq, respectively, 
the outer limb of the disc is about 1.1 Rq far from the 
secondary surface. Therefore, our results are realistic. 
According to the disc dimensions, it is a thick disc which 
its temperature is more than the temperature of the 
secondary and slightly less than the temperature of the 
primary. 

Altogether, on the basis of our results, AV Del has 
an accretion disc, while Mader et al. (2005) rejected ex- 
istence of activity in the system since they didn't found 
obvious signs of activity displayed by the system. How- 
ever, since there are so many free parameters of the disc, 
we can not claim that our fit is a unique solution (see 
Budaj et al. 2005), but existence of the disc in the sys- 
tem is certain. So, the period changes that were saw 
in the system (Qian, 2002) can be related to the mass 
transfer between two components which has created the 
accretion disc. 

Finally, Popper (1996) has noticed "although the two 
components appear to have comparable luminosities, 
the hotter and smaller component has the larger mass, 
with a mass ratio very provisionally about 1.9". So, 
according to the values we derived, our results confirm 
his results. 

According to our results, AV Del is a semi-detached 
system which the primary and secondary components 
filling are almost 49±1 and 100 percent of their respec- 
tive critical Roche lobes. 

Also, the absolute elements of both components de- 
termine the evolutionary state of the system. Enter- 
ing the results in the mass- luminosity (M-L), M^-mass 
(M-M), mass-radius (M-R) and H-R diagrams (the di- 
agrams are displayed in Fig. 6), it appears that AV 
Del contains a main-sequence primary star, while its 
secondary is on its way to the giant stage. 

Acknowledgements We wish to very thank Dr. 
Jano Budaj for his useful comments at various stage 
of the work. 
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• observed data 

theoretical light curve with disc 
theoretical light curve without disc 




Phase 



Fig. 4 The observed and the theoretical light curves of AV Del, while the theoretical light curves are drawn on the basis 
of the values derived via SHELLSPEC code. Solid circles show the observed data and the theoretical light curves are shown 
by continuous (with disc) and dashed (without disc) lines. 




Fig. 5 The Configuration of the binary system AV Del on the basis of the values derived via SHELLSPEC code. As we 
see, the primary star is surrounded by an accretion disc. 
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Fig. 6 M-R, H-R, M-L and M-M diagrams for the eclipsing binary system AV Del. The location of components are shown 
with star sign. 
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